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ABSTRACT 

We present a time, spectral and imaging analysis of the X-ray reflector in NGC 
4945, which reveals its geometrical and physical structure with unprecedented de- 
tail. NGC 4945 hosts one of the brightest AGN in the sky above 10 keV, but it is only 
visible through its reflected/scattered emission below 10 keV, due to absorption by a 
column density of ~4xl0 24 cm -2 . A new Suzaku campaign of 5 observations span- 
ning ^6 months, together with past XMM-Newton and Chandra observations, show a 
remarkable constancy (within <10%) of the reflected component. Instead, Swift-BAT 
reveals strong intrinsic variability on time scales longer than one year. Modeling the 
circumnuclear g clS clS cL thin cylinder with the axis on the plane of the sky, we show 
that the reflector is at a distance 30-50 pc, well within the imaging capabilities of 
Chandra at the distance of NGC 4945 (1"~18 pc). Accordingly, the Chandra imaging 
reveals a resolved, flattened, ^150 pc-long clumpy structure, whose spectrum is fully 
due to cold reflection of the primary AGN emission. The dumpiness may explain the 
small covering factor derived from the spectral and variability properties. 
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1 INTRODUCTION 

NGC 4945 is an almost edge-on (inclination angle ~ 80°) 
spiral galaxy hosting one of the nearest AGN (D ~3.7 Mpc, 
1"=18 pc, Mauersberger et al. 1996). NGC 4945 has been 
widely studied in the past, both in the soft and in the hard 
X-rays. It is the brightest Seyfert 2 and the second bright- 
est radio-quiet AGN after NGC 4151 in the 100 keV sky 
(Done et al. 1996), with a strongly absorbed (N H ~4x 10 24 
cm -2 , Itoh et al. 2008) intrinsic nuclear continuum, known 
to be extremely variable (Guainazzi et al. 2000, Madejski et 
al. 2000). Such a high column density completely blocks the 
primary nuclear emission below 8-10 keV while the emis- 
sion at higher energy is still visible, though heavily affected 
by Compton scattering and photoelectric absorption. A de- 
tailed analysis of the high energy variability (Madejski et 
al. 2000, Done et al. 2003, hereafter D03) already provided 
important (and, so far, almost unique) constraints on the 
geometry of the reflector: the variability above 10 keV im- 
plies that most of the observed emission is due to the pri- 
mary emission, and not to Compton scattering from other 
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directions (which would dilute the intrinsic variability, pro- 
ducing an almost constant observed emission). This, in turn, 
means that the reflector covers a solid angle < 10 deg as seen 
from the source (assuming a toroidal structure). This is also 
in agreement with the unusually low ratio between the ob- 
served, reflected emission below 10 keV, and the intrinsic 
flux in the same band as estimated from the emission above 
10 keV (/ <0.1%). Finally, a Chandra observation revealed 
a spatially resolved component, on a ~100 pc scale, whose 
spectrum is typical of Compton reflection from a ~neutral 
medium (D03). 

Here we present the analysis of new and archival Suzaku, 
XMM-Newton and Chandra observations, focusing specifi- 
cally on the determination, with unprecedented detail, of the 
dimensions and the geometrical structure of the circumnu- 
clear absorber/reflector in NGC 4945. A complete study, dis- 
cussing both the Chandra imaging data, and our new Suzaku 
observations, will be presented in a separate paper (Marin- 
ucci et al. 2012, in prep.). The work presented here expands 
the results summarized above, through two completely new 
studies: 

- the comparison between the intrinsic variability above 
10 keV (from Swift/BAT observations), and the reflection 
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variability below 10 keV (from a set of 8 XMM-Newton and 
Suzaku observations over a period of ~10 years); 
- a detailed analysis of the spatial extension of the reflec- 
tor, taking advantage of the full set of available Chandra 
observations, as discussed in the next Sections. 



2 OBSERVATIONS, AND DATA REDUCTION 

Chandra. NGC 4945 was observed by Chandra on 2000, 
January 27-28 for a total exposure time of 49 ks, with the 
AC1S camera. It was observed again four years later twice, 
for a total exposure time of 180 ks taking advantage of the 
HETG instrument. Data were reduced with the CIAO 4.3 
package (Fruscione et al. 2006) and using the Chandra Cal- 
ibration Data Base (CALDB) version 4.4.6, adopting stan- 
dard procedures. The imaging analysis was performed on 
a file including the two merged HETG-zeroth order im- 
ages reprojected on the ACIS-S image, applying the SER 
and smoothing procedures widely discussed in the literature 
(Tsunemi et al. 2001, Li et al. 2004, Wang et al. 2011), we 
therefore use a pixel size of 0.246 arcsec in Fig. [5] and 0.123 
arcsec in Fig. [4] After cleaning for background flaring we get 
a total of 34 ks and 173 ks for the 2000 and the 2004 HETG- 
zeroth order merged observations, respectively. Spectra were 
extracted on both the data set from three different regions: 
a circle with 25" radius (matching the XMM-Newton extrac- 
tion region, see below); a 12" x 6" box, and a central circle 
with 1.5 arcsec radius. We grouped all extracted spectra to 
have at least 20 total counts per new bin. 
Suzaku. The six Suzaku observations analysed in this paper 
were performed with the X-ray Imaging Spectrometer (XIS). 
The first one has been performed in 2006, has a duration of 
~100 ks, and has been presented by Itoh et al. (2008). The 
remaining five observations, with a duration of 40 ks each, 
have been performed between July 2010 and January 2011, 
as part of a monitoring campaign of this source. 

The event files were processed and calibrated adopt- 
ing standard procedures, discussed in several previous pa- 
pers (e.g. Maiolino et al. 2010). The source extraction ra- 
dius is 1.85' for all the 6 observations. Background spectra 
have been extracted from source-free regions with 3' radii. 
The 0.5-10 keV spectra extracted from the front-illuminated 
XIS0 and XIS3 have been co-added. 

XMM-Newton. We used 2 XMM- Newton observations 
performed with the EPIC-PN camera operated in large win- 
dow and medium filter modes. Source data 'cleaning' (exclu- 
sion of flaring particle background intervals) were performed 
with SAS 10.0.0 (Gabriel et al. 2004) via an iterative pro- 
cess which leads to a maximization of the Signal-to-Noise 
Ratio (SNR), similarly to that described in Piconcelli et 
al. 2004. Source spectra were extracted from a 25" radius 
circle; background spectra were extracted from source-free 
circular regions of the source field. 

Swift. The Swift/BAT 15-195 keV light curve has been ob- 
tained from the on-line catalog which at the time of our 
analysis was available for 65 months. 
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3 THE REFLECTOR GEOMETRY FROM 
TIME AND SPECTRAL VARIABILITY 

NGC 4945 offers a unique possibility to perform a compari- 
son between the variability of the intrinsic and the reflected 
X-ray emission, because it is the only know AGN with both 
(a) Compton thickness to>1, and (b) a strong intrinsic vari- 
ability above 10 keV. For all the other known bright AGN, 
such a detailed direct comparison between the variability of 
the two components is impossible because one of the two 
above mentioned conditions is not satisfied: 

- (a) if the AGN is not Compton-thick, the Compton reflec- 
tion component is diluted by the dominant direct compo- 
nent. The separation of the two components is never precise, 
and typical uncertainties in the flux of the reflected compo- 
nent are of the order of 20-30% even for long observations 
of the brightest AGN. 

- (b) the AGN is Compton-thick, but the intrinsic emission is 
either not visible (due to N H >10 25 cm" 2 , as in NGC 1068), 
or, anyway, not variable, probably because of the dilution 
due to multiple Compton scatterings (such as in the Circi- 
nus galaxy). 

NGC 4945 shows one of the brightest reflection dominated 
spectra below 10 keV, with flux measurement uncertainties 
of the order of ~1%, and a high-quality light curve of its in- 
trinsic emission above 10 keV, provided by Swift/BAT. The 
highest quality spectra below 10 keV have been obtained in 
the eight XMM-Newton and Suzaku observations presented 
in the previous Section. 

The 3-10 keV spectra of these observations have been 
fitted with a model consisting of a reflection continuum 
(PEXRAV, Madgziarz & Zdziarski 1995), three emission 
lines for neutral, He-like and H-like iron, and an additional 
power law, requested by the fit, and associated to the dif- 
fuse emission (see below). From the analysis of the HETG 
grating spectra we rule out the presence of emission features 
associated to the Iron K/3 line. The analysis has been per- 
formed with the XSPEC 12.7 code (Arnaud et al. 1996). 
All quoted errors are at the 90% confidence level for one 
interesting parameter. This model successfully reproduces 
all the individual spectra, with best fit parameters typical 
of reflection-dominated sources. The details of the spectral 
model have been discussed in Itoh et al. 2008, Schurch et 
al. 2002 and D03, and will be further analyzed in a forth- 
coming paper, together with the soft emission. What is in- 
teresting in this context is instead the possible variability 
of the reflected component. In this respect, the following re- 
sults are relevant: 

1) The analysis of the Chandra spectrum of the resolved 
12x6 arcsecond box, shows that the central emission is well 
fitted by a pure reflection continuum, with no need for an 
additional power law, so confirming that this component is 
needed only to reproduce the circumnuclear emission. Sim- 
ilarly, the XMM-Newton spectrum, extracted from the cen- 
tral 25", requires only a minor contribution from the addi- 
tional power law (about 1/10 of that in the Suzaku spectra). 

2) If we fit the eight spectra with constant reflection com- 
ponent and emission lines, leaving the additional power law 
as the only variable component, we obtain an equally good 
fit, with no significant residuals in any individual spectrum. 
This implies that formally, all the observed variability may 
be due to galactic sources. 
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Figure 1. Upper panel: Swift/BAT light curve of the 15-195 keV 
emission of NGC 4945, grouped in time bins of 1 month (black) 
and five months (red). Lower panel: expected maximum variabil- 
ity as a response of a cylindrical reflector to the BAT light curve, 
as a function of the distance form the central source. The value 
corresponding to the upper limit to the observed reflection vari- 
ability is ~35 pc. 



3) The whole luminosity within the central 1.85 arcmin re- 
gion is only ~2xl0 39 erg s _1 . It is therefore not surprising 
if significant variability due to single sources is observed. 
Based on the above considerations, we conclude that the 
best way to estimate the reflection variability is to fit the 
8 spectra with the model described above, with only three 
free parameters: the slope and the normalization of the addi- 
tional power law, and the flux F of the reflection (continuum 
plus lines) component, normalized to the value of the first 
interval, for ease of comparison. The results are reported in 
Table 1, and show no significant variability in any observa- 
tion. If we assume a constant value of the relative flux F, we 
obtain <F>=0.98±0.02, and a marginally better fit (the \ 2 
increases by 3, with 7 more degrees of freedom). The disper- 
sion around the central value is significantly smaller than 
the uncertainties on the individual values. This is due to the 
partial degeneracy among the spectral parameters. Based 
on these results, we estimate an upper limit to the observed 
reflection variability of ~4%, corresponding to twice the dis- 
persion, or, equivalently, to the ratio between the maximum 
and minimum value of F within the measured dispersion. 
This result can be compared with the variability observed 



Figure 2. Two-color (red:0.3-2 keV; green: 2-10 keV) Chandra 
image of the l'xl' central region of NGC 4945. 

at higher energies. The average 15-150 keV flux of NGC 4945 
is of the order of 5-10 mcrab, therefore the Swift/BAT mon- 
itoring provides reliable, high S/N light curves with reason- 
ably short time bins. In Fig. [1] we show the 65-month light 
curve, in bins of one and five months. In order to obtain a 
semi-quantitative estimate of the reflector distance, we mod- 
eled the circumnuclear gas as a cylinder with the axis on 
the plane of the sky, with radius R and height H, such that 
H/R<0.1, as implied by the analysis of D03. We then used 
the observed BAT light curve as an input, and determined 
the maximum observable flux variation of the reflected com- 
ponent as a function of the distance R, by adding the contri- 
bution of each element of the reflector. Qualitatively, we ex- 
pect that the more distant the reflector, the lower variability 
is observed in the reflected light curve, because the intrinsic 
variability is smoothed out by the different lengths of the 
light paths corresponding to each reflecting element. Quan- 
titatively, the result is shown in Fig. [Tp. In order to observe 
a variation not larger than 4% (reminding that <F>=0.96- 
1.00), the reflector must be at a distance larger than 35 pc. 
This result is obviously approximate, due to the assump- 
tions on the geometry, and to the limited sampling of the 
reflection fluxes, which may have caught the source in the 
same reflected states by chance. However, we believe it pro- 
vides a solid order of magnitude estimate of the distance of 
the reflector. We conclude noting that, considering that the 
Chandra resolution is ~ 0.25 arcsec, (~5 pc at the distance 
of NGC 4945) , the X-ray reflector should be easily resolved 
by Chandra, consistent with the conclusions of D03. 



4 THE REFLECTOR GEOMETRY FROM 
DIRECT IMAGING 

In Fig.[2]the l'xl' central region of NGC 4945 is shown. Af- 
ter running the script CHANDRA_repro on the three Chan- 
dra event files we merged them all with the merge_ALL 
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Figure 3. Spectra (rescalcd to a common effective area), model 
and A\ 2 from the observations. The A\ 2 values are shown for 
all the spectra in the central panel, while only some of the spec- 
tra are shown, for clarity. Upper panel: first Suzaku observation 
in 2010 (black); Suzaku in 2011 (red), Suzaku in 2006 (green); 
XMM-Newton in 2004 (blue). Lower panel: the first 2010 Suzaku 
observation is shown again (black) for ease of comparison; the re- 
maining spectra are from the Chandra 2004 data, extracted from 
a 25" radius (blue) a 12x6" box (red) and a 1.5" radius circle 
(green). It is worth noting that the spectra between 6 and 7 keV, 
dominated by the iron line emission, are almost identical in all 
observations, except for the last one from Chandra, and that the 
total continuum decreases with the area of the extraction region. 
Both these features are evidence that the variability is due to the 
diffuse galactic component, while the reflection remains constant, 
as quantitatively demonstrated by the spectral fits (Table 1). 



script and filtered them in energy between 0.3-2 keV (in red) 
and 2-10 keV (in green). A soft X-ray emitter can be seen 
in the top right region of the image, as already discussed 
in Schurch et al. 2002, where they called it the 'plume', 
reproduced by a single temperature thermal model (kT~ 
0.6 keV), interpreted as a mass-loaded superwind emanated 
from the central nuclear starburst. With the inclusion of the 
two 2004 observations in the image, we notice a further soft 
lobe in the bottom-left direction indicating a second, pos- 
sible emitter partially absorbed by a dust-lane crossing the 



nucleus, aligned with the clumpy structure shown in Fig. [4] 
For this analysis we refer to a future work (Marinucci et al., 
in prep.). 

The nuclear 2-10 keV emission is embedded in a 12" x 6" box, 
corresponding to a 180 pc x 90 pc region, which is consistent 
with the starbust ring traced by molecular gas (Moorwood 
et al. 1996, Marconi et al. 2000, Curran et al. 2001). 

We performed a spectral analysis of three zones, corre- 
sponding to the XMM-Newton extraction region, a 12" x 6" 
box, including most of the hard emission (Fig. [2]), and a cen- 
tral circular region with radius 1.5". The spectra, shown in 
Fig. [3] are fitted with the same reflection model as described 
above, with a contribution of the additional power law de- 
creasing from the larger to the small region. In particular: 
(a) the best fit model of the spectrum from the 25" circle is 
fully consistent (in each individual parameter) with that of 
the XMM-Newton observations; (b) the spectrum from the 
12" x 6" box is again fitted with the same model, but with 
a relative flux of the reflection component F=0.86±0.07, in- 
dicating that most of the reflection comes from this region, 
with only a remaining ~10% from outer regions; (c) the 
model from the central annulus is fitted with a pure re- 
flection component, with no need for the additional power 
law. The relative flux is F=0.58±0.06. This shows that the 
reflection emission is moderately concentrated in the pro- 
jected central ~ 30 pc, but with a significant contribution 
from the outer regions. In our simple cylindrical scheme, if a 
small (a few degrees) inclination angle is assumed, this con- 
centration in the central unresolved region is expected due 
to the projection effect, and is not related to a small physi- 
cal distance from the X-ray source. A direct, complementary 
way to demonstrate that the hard X-ray emission is domi- 
nated by the reflected, diffuse component is by comparing 
the image in the full hard X-ray spectral interval (2-10 keV), 
with the exclusion of the 6.2-6.7 energy range, with the "iron 
line image" obtained filtering only the 6.2-6.7 keV spectral 
interval (Fig. [4j) , which is clearly dominated by the iron emis- 
sion line (Fig. [3}- We then overimposed the contours of the 
hard X-rays on the Iron Ka emitter. The striking match be- 
tween these two images demonstrates that the cold reflection 
and the iron emission lines originate exactly from the same 
circumnuclear material and we are actually able to take a 
screenshot of the inner reflecting structure of NGC 4945. 

Finally, we notice that the emission in Fig. [4] originates 
from a non homogeneous structure, with visible clumps and 
empty regions with sizes of the order of tens of parsecs. The 
differences in counts among these sub-regions are highly sig- 
nificant (same area regions of 10x10 pixels contain from 
< 10 to >100 counts). A smaller-scale dumpiness of these 
structures, not resolved by Chandra, may explain the low 
covering factor (<10%) of the reflector, as estimated from 
the high intrinsic variability above 10 keV, and from the low 
ratio between the reflected and intrinsic components (Made- 
jski et al. 2000). 



5 CONCLUSIONS 

The time, spectral and imaging analysis of the X-ray 
emission of NGC 4945 based on XMM-Newton, Suzaku, 
Swift/BAT and Chandra data, characterizes of the AGN 
circumnuclear reflector with unprecedented detail. 
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Figure 4. Image of the central region in the line emission band (6.2-6.7 keV, left) and in the hard X-ray band (right). 



Par am. 


Value 


OBS 


Date 


F e 


r 


1.6±0.1 


Suzaku 1 


01/06 


1.00±0.06 


E(l) a 


6.39±0.01 


Suzaku 2 


07/10 


1.00±0.08 


EW(l) b 


1020±100 


Suzaku 3 


07/10 


0.99±0.08 


E(2) a 


6.66±0.01 


Suzaku 4 


07/10 


0.97±0.08 


EW(2) 6 


130±20 


Suzaku 5 


08/10 


0.98±0.08 


E(3) a 


7.02±0.01 


Suzaku 6 


01/11 


0.94±0.08 


EW(3)' 


140+ 20 

lw -10 


XMM 1 


01/00 


0.97±0.05 


N c 


3.7±0.7 


XMM 2 


01/04 


0.98±0.06 


r 2-10 


i.3xicr 12 


Chandra 25" 


05/04 


1.08±0.11 


T d 

^2-10 


2.2xl0 :i9 


Chandra BOX 


05/04 


0.86±0.07 


xVd.o.f. 


1308/1205 


Chandra 1.5" 


05/04 


0.58±0.06 



Table 1. Best fit values for the reflection component. a : emission 
line energies, in keV; : line equivalent widths, in eV; c : normal- 
ization in units of 10 — 3 ph cm" - 2 s" 1 keV" 1 ; d : 2-10 keV flux and 
luminosity of the reflection spectrum, in c.g.s. units. e : relative 
flux. 



1) The comparison between the strong intrinsic variability 
measured by Swift/ B AT and the constant reflection spectra 
from XMM-Newton and Suzaku observations over more 
than ten years implies a distance of the reflector D>35 pc. 

2) The Chandra image, obtained combining two observa- 
tions performed in 2000 and 2004, reveals an extended 
hard emission on projected scales of ~200xl00 pc. The 
spectrum of this emission is entirely reproduced by a cold 
reflection model. The central 30 pc account for about 50% 
of the whole emission. 

3) We show the first X-ray image of the inner reflector of 
an AGN. The large scale structure of the emission region is 
clumpy and asymmetric, and the central region is resolved 
by Chandra. A smaller-scale, unresolved dumpiness may 
explain the low covering factor of the reflector inferred by 
the high variability of the intrinsic emission at E>10 keV, 
and by the low ratio between reflected and intrinsic spec- 
trum (Madejski et al. 2000, D03). 
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